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' Abstract 

We study the classical and quantum models of a scalar field Friedmann-Robertson- Walker 
£j " (FRW) cosmology with an eye to the issue of time problem in quantum cosmology. We intro- 

duce a canonical transformation on the scalar field sector of the action such that the momentum 
conjugate to the new canonical variable appears linearly in the transformed Hamiltonian. Using 
this canonical transformation, we show that, it may lead to the identification of a time parameter 
for the corresponding dynamical system. In the cases of flat, closed and open FRW universes the 
classical cosmological solutions are obtained in terms of the introduced time parameter. Moreover, 
this formalism gives rise to a Schrodinger- Wheeler-DeWitt equation for the quantum-mechanical 
description of the model under consideration, the eigenfunctions of which can be used to construct 
the wave function of the universe. We use the resulting wave functions in order to investigate the 
possible corrections to the classical cosmologies due to quantum effects by means of the many- 
worlds and ontological interpretation of quantum cosmology. 
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00 ■ 1 Introduction 

o: 

Since 1967 when canonical quantum theory of gravity was first introduced by DeWitt [1], many efforts 
have been made in this area and the corresponding results have been followed by a number of works, 
the main motivations of which lie in the results of the unification theories and cosmology [2] . The main 
object of such a theory, i.e., the quantum state of the gravitational system, is then described by a wave 
function, a functional (on superspace) of the 3-geometry and matter fields presented in the theory, 
satisfying the Wheeler-DeWitt (WDW) equation [3]. One of the most important features related to 
the WDW wave function is a problem which is known as the "problem of time". Indeed, unlike the 
case of usual quantum theory, the wave function in quantum gravity is independent of time, which is 
a reflect to the fact that general relativity is a parameterized theory in the sense that its (Einstein- 
Hilbert) action is invariant under time reparameterization. Because the existence of this problem, 
the canonical formulation of general relativity leads to a constrained system, and its Hamiltonian is a 
superposition of some constraints, the so-called Hamiltonian and momentum constraints. A possible 
way to overcome this problem is that one first solves the equation of constraint to obtain a set of 
genuine canonical variables with which one can construct a reduced Hamiltonian. In this kind of time 
reparameterization, the equations of motion are obtained from the reduced physical Hamiltonian and 
describe the evolution of the system with respect to the selected time parameter [4]. 

The problem of time was first addressed in [1] by DeWitt himself. However, he argued that the 
problem of time should not be considered as a hindrance in the sense that the theory itself must 
include a suitable well-defined time in terms of its geometry or matter fields. In this scheme time is 
identified with one of the characters of the geometry, usually the scale factors (in cosmological models) 
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of the geometry and is referred to as the intrinsic time, or with the momenta conjugate to the scale 
factors, or even with a scalar character of matter fields coupled to gravity in any specific model, 
known as the extrinsic time. Identification of time with one of the dynamical variables depends on 
the method we use to deal with theses constraints. Different approaches arising from these methods 
have been investigated in detail in [5] and [6]. 

In this paper, we first deal with a FRW cosmology with a scalar field minimally coupled to the 
gravity. The phase-space variables of such a model turn out to correspond to the scale factor a(t) of the 
cosmological model and a scalar field (p(t) with which the action of the model is augmented. We then 
introduce a canonical transformation on the scalar field sector of the action such that in terms of the 
new canonical variables, the transformed Hamiltonian contains a linear momentum. This procedure 
causes the Hamiltonian takes the form of a Schrodinger one in which the variable its conjugate 
momentum appears linearly in the Hamiltonian plays the role of a time parameter. Moreover, this 
formalism gives rise to a Schrodinger- Wheeler-DeWitt (SWD) equation for the quantum-mechanical 
description of the model under consideration, the eigenfunctions of which can be used to construct 
the wave function of the universe. In the cases of flat, closed and open FRW models we present the 
classical and quantum cosmological solutions in terms of the introduced time variable and investigate 
the behavior of the corresponding universe in each case separately. 



2 The classical model 



In this section we consider a FRW cosmology with a scalar field with which the action of the model 
is augmented. In a quasi-spherical polar coordinate the geometry of such a space-time is described 
by the metric 

l-kr 2 +7 " 2 (^ 2+sin2 ^) ' (!) 



ds z 



-N 2 (t)dt 2 + a 2 (t) 



where N(t) is the lapse function, a(t) the scale factor and k = 1, and —1 corresponds to the closed, 
flat and open universe respectively. Since our goal is to study a procedure in which how a variable 
may play the role of a time parameter, we do not include any matter contribution in the action. Let 
us start from the action (we work in units where c = h = 16nG = 1) 



/ 



d XyJ—g 



R - ^d^d^ - v(<f>) 



(2) 



where g is the determinant of the metric tensor g^, R is the Ricci scalar corresponding to g^ v and 
V(4>) is the potential function for the scalar field <f)(t). By substituting (1) into (2) and integration 
over spatial dimensions, we are led to a point-like Lagrangian in the minisuperspace {a, <f\ as 

C = -3ad 2 + ^a 3 <p 2 + 3ka- a 3 V((f>), (3) 

in which we have set N = 1 so that the time parameter t is the usual cosmic time. The momenta 
conjugate to each of the above variables can be obtained from their definition as 



P - d A 

oa 



-6aa, 



dC 31 

d<t> 



(4) 



In terms of these conjugate momenta the canonical Hamiltonian, which is constrained to vanish, is 
given by 



H 



Now, we consider the following canonical transformation 
of the action [7] 



, Pa, 



p2 

2 ' 



(5) 

(T,Pr) on the scalar field sector 

(6) 
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Under this transformation Hamiltonian (5) takes the form 



n 



^- + ^-3ka + a 3 V(T). 
12a a 6 



(7) 



We see that the momentum Pt appears linearly in the Hamiltonian. This means that the parameter 
T may be interpreted as a time parameter when one deals with the quantum version of the model. 
Considering the parameter T as a clock parameter, the Hamiltonian (7) will be a time-dependent 
function. Such Hamiltonian describes a system which exchanges energy with the surrounding en- 
vironment. However, in the case of cosmology where the system under consideration is the whole 
universe, a surrounding environment does not have any meaningful interpretation. Therefore, such a 
Hamiltonian and the corresponding time parameter do not seem to be suitable unless the potential 
function is constant which for the sake of simplicity we consider it to be zero. In this case the classical 
dynamics is governed by the Hamiltonian equations, that is 

a = { a ,H} = -^ 



6a ' 



Pa = {Pa,-H} = - 

T = {T,H} = \, 



+ +3k, 

Via 2 - a 4 ' 



(8) 



{ P T = {P T ,H} = 0. 

Now, we would like to rewrite the above equations in terms of the new time parameter T. By definition 
l' = Wr an d w ith the help of the third equation of the above system, the classical equations of motion 
can be rewritten as follows 



a' = -\a 2 P a , 

P' a = -^P 2 a + Z -t + ^a\ 



(9) 



where we take Pt = A = const, from the last equation of (8). We also have the constraint equation 
% = from which we obtain P 2 = — 36/ca 2 . Eliminating P a between this equation and the first 
equation of (9) results 

a 12 = AV - ka 6 , (10) 

where A 2 = A/3. In the following we present the solutions to this equation according to the various 
values of the curvature index k. 

• The flat universe: k = 0. In this case equation (10) admits the solutions 



a(T) = aoe 



(11) 



where ao is an integration constant and the positive (negative) sign corresponds to an expanding 
(contracting) universe. For a positive sign in the power of the exponential term, the evolutionary 
behavior of the corresponding universe based on (11) is like a de Sitter universe, i.e., begins with a 
zero size at t = — oo and follows the exponential law expansion at late time of cosmic evolution. For 
a negative sign, on the other hand, the behavior is opposite. The universe decreases its size from 
large values of scale factor at t = — oo and tends to zero at the late time. In summary, what we have 
shown is that a flat FRW model in the above gauge for the time parameter has two separate classical 
solutions, one describes an expanding universe while another represents a contracting one. In figure 
2 (the left figure) these two classical scale factors are plotted for typical values of the parameters. 
• The closed universe: k = 1. In this case equation (10) takes the form 



,'2 



AV 



(12) 
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It is clear from this equation that the scale factor is restricted to belong to the interval < a < VX 
Performing the integration we get the following implicit relation between T and a(T) 

a 2 (T) 

_ a„±2\T 



A 2 + Av/A 2 -a 4 (T) 



Ae ±2X1 , (13) 



where A is an integration constant. It is seen that the above mentioned restriction on the values of 
the scale factor shows itself in the expression under square root. Solving (13) for a(T), we get 

V 1 Vl + X 2 A 2 e ±iXT V ' 

Again, we see that the classical closed model has also two separate solutions (one for the positive sign 
in the power of the exponential function and another for the negative sign) which we have plotted 
them in figure 4 (the left figure). As is clear from the figure, in both of these models a contraction 
period is followed by an expansion era. However, the turning point in which the transition between 
two expanding and contracting phases is occurred, is not the same for them. 
• The open universe: k = — 1. In this case we have 



a' 2 = AV + a 6 , (15) 



where upon integration we obtain 



_ a„±2\T 



A 2 + AVA 2 + a 4 (T) 

from which we have 



Ae ±2M , (16) 



,(T) ^ ±XT (17) 
V ' Vl - X 2 A 2 e ±iXT 

These solutions are plotted in the left figure of figure 5. As the figure shows, one branch of the 
classical solutions begins with zero size at t = — oo, increase the size as time grows and tends to a 
singularity at which the scale factor blows up. The another branch has the opposite behavior, i.e., 
begins with a singularity at which the scale factor has a large size and then follows a contracting 
epoch and finally reaches the zero size as t — > +oo. 



3 Quantization of the model 

In this section we would like to see that how the classical picture will be modified when one takes 
into account the quantum mechanical considerations in the problem at hand. Before going to do 
this, a remark about the canonical transformation which led us to the Hamiltonian (7) is in order. 
The canonical transformation (6) is applied to the classical Hamiltonian (5), resulting in Hamiltonian 
(7) which we are going to quantize. To make this acceptable, one should show that in the quantum 
theory the two Hamiltonians are connected by some unitary transformation, i.e. the transformation 
(6) is also a quantum canonical transformation. A quantum canonical transformation is defined 
as a change of the phase-space variables (q,p) — > (q',p') which preserves the Dirac bracket, i.e. 
[QiP] = i = [q'(q,p),p'(q,p)]- Such a transformation is implemented by a unitary operator C such 
that q'(q,p) = CqC^ 1 and p'(q,p) = CpC^ 1 [8]. Under the act of this unitary transformation state 
vectors will be transformed as tp' = Cip and any physical observable A(q,p), being an operator on 
the Hilbert space of states, transforms according to A' = CAC~ X . Now, it is easy to see that if C 
is a unitary operator (C* = C _1 ) associated with some quantum canonical transformation, then the 
Hermitian property does not change after it transforms things from one coordinate system to the 
other. Indeed, let A = A^ be a Hermitian operator associated to some physical observable. Under the 
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act of a canonical transformation C, we have A'^ = (CAC^ 1 )^ = (C^A^C* = CAC' 1 = A', which 
means that A' is also a Hermitian operator 1 . For our case it is easy to see that the canonical relations 
[<f>, P,/)] = i yield [T,Pr] = 4>P^ X ^\P^ = i. This means that the transformation (6) preserves the 
Dirac brackets and thus is a quantum canonical transformation. Therefore, use of the transformed 
Hamiltonian (7) for quantization of the model is quite reasonable. 

Now, our starting point to quantize the model is to construct the WDW equation fi^(a,T) = 0, 
in which % is the Hamiltonian operator where its classical expression is given by (7) and \I/(a,T) 
is the wave function of the universe. With the replacement P a — > —i-§^ and similarly for Pt and 
choosing the potential to be zero, the WDW equation reads 

(02 d d \ 

a2 d^ + m d~a ~ 121 df ~ 3Qka ) * (a ' T) = °' (18) 

where the parameter s represents the ambiguity in the ordering of factors a and P a in the first term 
of (7). This equation takes the form of a Schrodinger equation and we separate its variables as 

*(o,T) = e iET ^(a), (19) 

leading to 

r _ # d .1 

V>(a) = 0, (20) 



d 2 d 
a 2 —, + sa— + 12(E - 3ka 4 ) 
da z da 



where E is a separation constant. In the following, as in the classical model, we shall deal with the 
solution to this equation in k = ±1,0 cases separately. 

• The flat universe: k = 0. In this case equation (20) has the following solutions 

$ E {p) = a sin (2V3E In a) + c 2 cos (2,V3E\n aj , (21) 

where c\ and C2 are the integration constants and we have chosen s = 1 . If, without losing the general 
character of the solutions, we set C2 = the eigenfunctions of the SWD equation can be written as 

^ E (a,T) = e iET sin(2v / 3^1na) . (22) 

We may now write the general solution to the SWD equation as a superposition of its eigenfunctions, 
that is 

*(o,T)=/ A(E)V E (a,T)dE, (23) 
Jo 

where A(E) is a suitable weight function to construct the wave packets. By using the equality 



[°° e^ x sin y/Wr = ^l e -(™/^) (24) 
Jo 2y' 2 

we can evaluate the integral over E in (24) and simple analytical expression for this integral is found 
if we choose the function A(E) to be a quasi-Gaussian weight factor A(E) = e~ lE (7 is an arbitrary 
positive constant), which results in 

*(o, T) = J™ e~~ tE e lET sin (2V3E In a) dE. (25) 



x At this step a remark is in order about the self-adjointness problem in the relation T = <j>/P$. In general, the 
self-adjointness of <f> and P$ does not lead to this property for T, and with a non-self-adjoint operator the construction 
of a self-adjoint Hamiltonian will be complicated. However, note that what is explicitly appeared in the T-sector of the 
Hamiltonian (7) (with zero potential term) is Pt (which is clearly a self-adjoint operator if P^ is) and not T. Also, in 
the a-sector of the Hamiltonian the self-adjointness is encoded in the chosen factor ordering, see equation (18) below. 
Therefore, we may claim that our theory is based on a self-adjoint Hamiltonian. Since with the Hamiltonian (7), the 
quantization of the model yields a Schrodinger-like equation, in what follows, the role of the variable T looks like the 
time parameter in Schrodinger equation in terms of which the dynamics of the wave function and other physical variables 
should be obtained. 
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Figure 1: The figures show the square of the wave function (26) and its contourplot. 



Using of the relation (24) leads to the following expression for the wave function 



^(a,T) = TV 



In a 



31n 2 a' 



( 7 - iTfl 2 



exp 



7 



iT 



(26) 



where J\f is a numerical factor. 

In figure 1 we have plotted the square of the wave function and its contourplot for typical numerical 
values of the parameters. As this figure shows, at T = 0, the wave function has a dominant peak in 
the vicinity of some non-zero value of a. This means that the wave function predicts the emergence 
of the universe from a state corresponding to its dominant peak. As time progresses, the wave packet 
begins to propagate in the a-direction, its width becoming wider and its peak moving with a group 
velocity towards the greater values of scale factor. The wave packet disperses as time passes, the 
minimum width being attained at T = 0. As in the case of the free particle in quantum mechanics, 
the more localized the initial state at T = 0, the more rapidly the wave packet disperses. Therefore, 
the quantum effects make themselves felt only for small enough T corresponding to small a, as 
expected and the wave function predicts that the universe will assume states with larger a in its 
late time evolution. Now, having the above expression for the wave function of the universe, we are 
going to deal with the question of the recovery of the classical solutions, in other words, how can the 
Wheeler-DeWitt wave functions predict a classical universe. In this approach, as we have done above, 
one usually constructs a coherent wave packet with good asymptotic behavior in the minisuperspace, 
peaking in the vicinity of the classical trajectory. On the other hand, to show the correlations between 
classical and quantum pattern, following the many- worlds interpretation of quantum mechanics [9], 
one may calculate the time dependence of the expectation value of a dynamical variable q as 



< q > (t) 



< ^\q\$! > 

< > 



(27) 



Following this approach, we may write the expectation value for the scale factor as 

J™^*(a,T)a^(a,T)da 



<a> (T) 



f™**(a,T)*(a,T)da 



which yields 



< a > (T) ~ 



T 2 + 7 2 + 3 7 



s 7 



(28) 



(29) 



T 2 + 7 2 + 12 7 

It is important to classify the nature of the quantum model as concerns the presence or absence of 
singularities. For the wave function (26), the expectation value (29) of scale factor never vanishes, 
showing that these states are nonsingular. Indeed, the expression (29) represents a bouncing universe 
with no singularity where its late time behavior coincides to the late time behavior of the classical 
solution. This means that the quantum structure which we have constructed has a good correlation 
with its classical counterpart. 
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The issue of the correlation between classical and quantum schemes may be addressed from another 
point of view. It is known that the results obtained by using the many-world interpretation agree 
with those that can be obtained using the ontological interpretation of quantum mechanics [10]. In 
Bohmian interpretation, the wave function is written as ^(a,T) = Q(a,T)e zS ( a ' T \ where f2(a, T) and 
S(a,T) are real functions, which using the expression (26), a simple algebra gives them as 

n(a, T)=M ( 7 2 + T 2 ) ~ 3/4 exp (-37 ^° 2 ) In a, (30) 

x 3 T 3Tln 2 a , 

S(a, T ) = 2 arctan ~ ~ 7 2 + T 2 • ( 31 ) 

In this interpretation the classical trajectories, which determine the behavior of the scale facto is 
given by P a = || , where by means of the relation (4) reads 

I da T 

--r= = — —, 7 In a, (32) 
adT -f 2 + T 2 y ' 

which, after integration we get 



a(T) = a e^ 2+T2 . (33) 

This solution has the same behavior as the expectation values computed in (29) and like that is free 
of singularity, figure 2 shows the behavior of the classical scale factor (11), the quantum mechanical 
expectation value of the scale factor (29) and its Bohmian version (33) versus time for some typical 
numerical values of the parameters. We see that instead of two separate classical solutions, the 
quantum models predict a bouncing universe in which the universe decreases its size, reaches a 
minimum and then expands forever. The origin of the singularity avoidance may be understood by 
the existence of the quantum potential which corrects the classical equations of motion. According to 
the Bohm-de Broglie interpretation of quantum mechanics and also its usage in quantum cosmology, 
upon using the polar form of the wave function ^(a, T) = Q.(a, T)e tS ^ a ' T \ in the corresponding wave 
equation, we arrive at the modified Hamilton-Jacobi equation as 

n(qi,Pi = ^-)+Q = Q, (34) 



where Pj are the momentum conjugate to the dynamical variables % and Q is the quantum potential. 
For the wave equation (18) this procedure leads 

1 fdS\ 2 . 1 fdS\ 



na\8a) + -Afff)- Zka + Q = ^ <35) 
in which the quantum potential is defined as 

_j_dn 

12afi da 2 12a 2 n da' [ ' 

Now, inserting the relation (33) in (36), we can find the quantum potential in terms of the scale factor 
as 

\ ao ' 

As this relation shows, the potential goes to zero for the large values of the scale factor. This 
behavior is expected, since in this regime the quantum effects can be neglected and the universe 
evolves classically. On the other hand, for the small values of the scale factor the potential takes a 
large magnitude and the quantum mechanical considerations come into the scenario. This is where the 
quantum potential can produce a huge repulsive force, which may be interpreted as being responsible 
of the avoidance of singularity. 
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Figure 2: From left to right: The classical scale factors (11), the expectation value (29) and the Bohmian value (33). 




Figure 3: The figures show the square of the wave function (43) and its contourplot. 



• The closed universe: k = 1. Equation (20) in this case admits the following solutions in terms 
of the modified Bessel functions 



ip E (a) = ciK i ^ E (3a 2 ) + c 2 I i ^ E (3a 2 ) . 



(38) 



We impose the boundary condition lim; C _ s . 00 ip(x) = on the wave function which restricts us to 
consider the modified Bessel function K as solution. Thus the eigenfunctions of the SWD equation 
take the form 



which their superposition gives the general solution as 

roo 2 



3a' 



$(a,T) = J™ -vC{y)eh iv * T K lv (3a 2 ) du, 



(39) 



(40) 



where v = V3E. Since the function Ki v {z) decreases when v grows, thus the small i/'s have dominate 
contribution to the above integral. This allows us to use the approximation v 2 /2> ~ v. Therefore, 
with the choice of the weight factor as C(y) = |e 7i/ we get 



*(o,T) 

Now, by using the equalities [11] 



Ki v (3a 2 ) dv. 



(41) 



cos(Pv)K iu (a)dv = Je-^ 0811 / 3 
2 



POO m- 

, / v wx{Bv)K iv {a)dv = -a sinh/3 e - acosh/3 , (42) 
Jo 2 



for a > and |Im/3| < ^, we can evaluate the integral in (41) to achieve the following analytic form 
for the wave function 



*(a, T)=N 1 + 3a 2 sin (7 + %T) 



-3a 2 cos(7+iT) 



(43) 



where again J\f is a numerical factor and 7 should be chosen as I7I < j. In figure 3 the wave 
function and its contourplot is plotted. This figure shows the universe expands and after reaches a 
maximum recollapses. It is clear that the two phases of cosmic evolution are connected to each other 



..: I;. 'I. n.!. m: 




Figure 4: From left to right: The classical scale factors (14), the expectation value (44) and the Bohmian trajectory 
obtained from (47). 



without meeting any kind of singularity. Now, with the help of the definition (27) we can evaluate 
the expectation value of the scale factor when the corresponding quantum universe is in the state 
(43). This calculation results 



< a > (T) = cosh -1 / 2 T x 
8 (e T + e 3T + e 5T ) + (l + e 2T ) 3 cos 7 + 4e T (l + e 4T ) cos 2 7 - (l + 3e 2T + 3e 4T + e 6T ) cos 37 
lle T + 16e 3T + lle 5T + (1 + e 2T ) 3 cos 7 + 2 (4 + 5e 2T + 4e 4T ) cos 2 7 - (1 + 3e 2T + 3e 4T + e 6T ) cos 3 7 



,(44) 



where is plotted in figure 4. Also, to extract the Bohmian trajectories of this model we may write 
the wave function (43) as \I/(a,T) = Q(a, T)e iS ^ a,T ^ in which 

n (a, T) = [l + 9a 4 (sin 2 7 cosh 2 T + cos 2 7 sinh 2 T) + 6a 2 sin 2 7 cosh 2 t] ^ e~ 3a2 cos 7 cosh T , (45) 
and 

2 m 3a 2 cos 7 sinh T 

b(a,T) = 6a sm^smhT + arctan ^ — -. (46) 

1 + 3a z sin 7 cosh T 

Substitution of the relation (46) in P a = ^ and using of the definition P a in (4), we obtain 



d a a 3 sinh T cos 7 + sin 7 — |a 4 (cos 27 — cosh 2T) sin 7 + 6a 2 cosh Tsin 2 7 



dT 1 — I a 4 (cos 27 — cosh 2T) + 6a 2 cosh T sin 7 



(47) 



This equation does not seem to have analytical solution. However, in figure 4, employing numerical 
methods, we have shown the qualitative behavior of a{T) for typical values of the parameters and 
initial conditions. As is clear from the figure, the scale factor repeats its expansion and contraction 
behavior as its expectation value also shown in this figure. We see that the quantization of the closed 
FRW cosmology predicts an evolution pattern for the corresponding universe in completely agreement 
with its classical dynamics which means that there is a good correlation between the classical and 
quantum models in this case as well. 

• The open universe: k = — 1. The solutions to the equation (20) with k = — 1 read as 

rl> E {a) = Cl J iV ^(3a 2 ) + c 2 y^(3a 2 ), (48) 

where J and Y are the Bessel functions. Removing the function Y from the solutions because of 
applying the same boundary condition as in the closed case, we are led to the following eigenfunctions 

y E (a,T) = e lET J iV ^(3a 2 ). (49) 

We may now write the general solutions to the SWD equation as a superposition of the above 
eigenfunctions 

^{a,T)= C{E)e iET J iV ^{3a 2 )dE. (50) 

J 

Unfortunately, unlike the previous two cases, it is not possible to find an analytical closed form for the 
wave function in the case of the open universe. Therefore, we should rely on the numerical methods 
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Figure 5: From left to right: The classical scale factors (17), the square of the wave function (50) and its contourplot. 

by using of the shifted Gaussian weig ht function C(E) = e~^ E -^ 2 . In fi gure 5 we have plotted 
the square of the wave function and its contourplot. The discussions on the comparison between 
quantum cosmological solutions and their corresponding form from the classical formalism are the 
same as previous models, namely the flat and closed models. Similar discussion as above would be 
applicable to this case as well. 

4 Summary 

In this letter we took a look at the problem of time in the framework of a scalar field cosmology. To 
do this, we have examined a canonical transformation (cj), Ps) — > (T,Pt) on the scalar field sector 
of the action which turns out to correspond to a specific choice of time parameter. In terms of this 
time parameter, we have solved the classical field equations in the cases of the flat, closed and open 
FRW universe and obtained the analytical expressions for the corresponding scale factors. As for the 
quantum version of these models, the use of the above mentioned canonical transformation allowed 
us to obtain a SWD equation in which the variable T plays the role of time parameter. In the cases 
of the flat and closed FRW models, we obtained exact solutions of the SWD equation and for the 
open model the wave function is represented numerically. The wave functions of the corresponding 
universes consist of some branches where each may be interpreted as part of the classical trajectories. 
We saw that since the peaks of the wave functions follow the classical trajectories, there seems to be 
good correlation between the corresponding classical and quantum cosmology. 

One of the main features of the quantum solutions in our presented model is that the singular 
region of the classical cosmologies is replaced by a bouncing period. In such a behavior which is 
shown in figures 1-5 the scale factor bounces (falls) from the contraction (expansion) to its expansion 
(contraction) eras. In addition to singularity avoidance, the appearance of bounce in the quantum 
model is also interesting in its nature due to prediction of a minimal size for the corresponding 
universe. We know the idea of existence of a minimal length in nature is supported by almost all 
candidates of quantum gravity. In this sense our results are qualitatively comparable with other 
works in which different quantization schemes are considered. We may name some of such works, for 
instance, as: 

• In [12], using the polymer quantization method of a k = FRW model, the volume of the 
corresponding universe is obtained as \/7 2 A 2 + 9i 2 , where 7 is the Barbero-Immirzi parameter and A 
is a parameter associated to the fundamental granularity of quantum geometry. 

• In [13], using a gauge-fixed Lagrangian to quantize a perfect fluid FRW cosmology, a bouncing 

scale factor (7 2 + i 2 ) 3 (" +1 ), for ordinary perfect fluid and a falling one [7 2 + (to — t) 2 ] 3 <"+ 1 ) for 
phantom perfect fluid is obtained, where 7 is an arbitrary positive constant and cj is the equation of 
state parameter of the perfect fluid. 

• In [14], a perfect fluid FRW cosmology is quantized with the use of the Schutz' representation 
of the perfect fluid which is shown may lead to an identification to a time parameter. The bouncing 
expression such as [^—^—t 2 + 1] 3 (!-") is obtained for the scale factor in this model. 



10 



• In [15] several bouncing solutions of various cosmological models are investigated and the mech- 
anisms behind the bounce are discussed. 
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